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We address the influence of realistic disorder on the effective magnetic monopole that is induced 
near the surface of an ideal topological insulator (Tl) by azimuthal currents which flow in response 
to a suddenly introduced external electric charge. We show that when the longitudinal conductivity 
o'xx ~ g{e^ /h) is accounted for, the apparent position of a magnetic monopole initially retreats from 
the TI surface at speed vm ~ oicg, where a is the fine structure constant and c is the speed of light. 
For the particular case of TI surface states described by a massive Dirac model, we further find that 
the temperature T = Hall currents vanish once the surface charge has been redistributed to screen 
the external potential. 
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Introduction- When a time-reversal-symmetry breaking 
perturbation opens a gap in the surface state spectrum of 
a three-dimensional topological insulator (TI) [U [5] , sur- 
face Hall currents and orbital magnetism are induced by 
electrical perturbations. This magneto-electric coupling 
effect can be attractively described [3] by adding a. E ■ B 
term to the electromagnetic Lagrangian. The duality of 
the resulting axion electrodynamics model |3] leads to a 
curious topological magneto- electric effect [SHZ] in which 
an electric charge placed above the TI surface induces 
Hall currents and associated orbital magnetization that 
appears to emanate from a magnetic monopole below the 
surface. 

In this paper we show that a non-zero TI surface state 
longitudinal conductivity a^x = an omnipresent 

experimental reality that is not captured by the axion 
electrodynamics model, qualitatively alters the topolog- 
ical magneto-electric effect. We find that when the ex- 
ternal charge is placed more than a screening length A 
from the surface, the monopole moves away with velocity 
Vm = acg. In the long-time limit the screened external 
potential becomes static. In this case we find that the 
orbital magnetization response depends on details of the 
surface state electronic structure, and that it vanishes 
in the particular case of a two-dimensional massive Dirac 
model with temperature T = and a Fermi level position 
outside the gap. 

Macroscopic Theory- We assume here that the TI surface 
has a well defined surface Hall conductivity and diffusion 
constant; this assumption can fail for very well developed 
quantum Hall effects. We first consider the limit in which 
the separation d between the external charge and the TI 
surface is larger than the screening length A. We intro- 
duce an external charge Qe located a distance d from the 
TI surface; since we wish to treat this object as a source 
of macroscopic inhomogeneity rather than as a contribu- 
tion to the disorder potential we imagine that Q ^ 1 and 
that d is longer than microscopic lengths. Currents flow 
in the TI surface in response to the electric fields from the 



external charge and the screening charges that accumu- 
late in the TI surface layer. Working in two-dimensional 
momentum space and assuming that the total electric 
field changes sufficiently slowly with time, we use the 
continuity equation to conclude that 

-^ = -27raxxqiQeM-qd)+nl^). (1) 

where d is the distance from the surface to the exter- 
nal positive charge Qe and n^^ is the Fourier transform 
of the induced surface state density. In Eq.([l]) we ne- 
glect the diffusion current, which is permissible at long 
distances as we show below. If we assume that the ex- 
ternal charge is introduced suddenly at time t = and 
that the two-dimensional (2D) density evolves in time in 
accordance with Eq. ([T]) we find that 

= -Qcxp{-qd){l - cxp{-qvMt)) (2) 

and that the total potential from external and screening 
charges is 

27re(5 

<Ptot{q,t) = exp{-q{d + VMt)). (3) 

q 

Here the monopole velocity vm = '^t^ctxx = c^cg is large 
unless the dissipative conductivity is much smaller than 
the quantum unit of conductance, i.e. unless the quan- 
tum Hall effect on the TI surface is very well developed. 
The potential at time t, which controls the instantaneous 
Hall currents and hence the instantaneous magnetization 
is identical to that from a external charge that is located 
not at vertical position d, but at vertical position d+VMt- 
As shown elsewhere|5j, because of the magneto-electric 
duality of axion electrodynamics, these Hall currents give 
rise to a magnetization that is identical to that produced 
by a magnetic monopole located at a distance d + VMt 
below the TI surface. Currents flow until macroscopic 
electric fields vanish. The topological magneto-electric 
effect is therefore purely transient in the d^ X limit. 
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Screening in the Quantum Hall Regime- This result can 
be extended by including the diffusion contribution to 
the surface current: 



Ik 



= -2TTa,j:xq{Qe 



— qd 



(4) 



The longitudinal conductivity, Uxx is related to the dif- 
fusion coefficient via the usual Einstein relation Uxx = 
[dn/dfije^Dp- Solving Eq. Q, we obtain the final ex- 
pression for the total electric potential on the surface: 



(l)tot{q,t) = 27reQe-«'^ 



1 + (gATF)"^e-('°^«'+27r^--«)*^ 



q + X 



TF 



(5) 

where Aj^]^ — 2'Kaxx/ Dp — 2TTVpe^ is the screening 
wavevector and Xtf the screening length. The longi- 
tudinal currents vanish for t — >■ oo due to the Einstein- 
relation cancellation between drift and diffusion contri- 
butions. The total potential for t — oo reduces to the 
standard result for Thomas- Fermi screening in 2D. Be- 
cause dn/d^L becomes extremely small when the quantum 
Hall effect is well developed, A can be much larger than 
typical microscopic length scales. 

Since the external potential remains large for i —> oo at 
length scales smaller than A, there will be a macroscopic 
orbital magnetic response to the screened potential if the 
contributions to the Hall current from the screened elec- 
tric field and from the induced density inhomogeneities 
do not cancel. Is there an Einstein relation for Hall cur- 
rents? Below we use a quantum kinetic theory to answer 
this question microscopically. We conclude that the an- 



swer is no in general. Both drift and diffusion type terms 
do appear. The contribution to the Hall current from 
density inhomogeneities can be understood as being due 
to non-uniform internal magnetic moment [S] densities. 
Moreover, for the two-dimensional massive Dirac equa- 
tion that is normally used to model TI surface states, 
the drift and diffusion Hall currents do cancel in the clean 
limit. 

Microscopic Theory- In the presence of an external po- 
tential the surface states of a 3D strong topological in- 
sulators can be described [TJ [2] approximately by a 2D 
massive Dirac Hamiltonian: 



MD 



rf^r^t (Bpfi + e0ext + C/dis) *. 



(6) 



Here Bp = [vpx,vpy, l^) is a p-dependent effective Zee- 
man field which acts on spinful surface electrons. With 
this choice for Bp, the Pauli matrices correspond to spins 
rotated by 7r/2 around the z axis, which we have taken to 
be normal to the surface. The mass term A breaks time- 
reversal-symmetry and is normally thought of as arising 
from proximity exchange coupling to an insulating ferro- 
magnet. For definiteness and without loss of generality, 
we take A > 0. C/dis describes an atomic scale disorder 

potential which we take to be created by short-range im- 
purities with concentration Uimp'- C/dis = ~ '"O- 
From now on we work in the system of units with h = 1. 

In order to address the transport properties of this 
model, we use a quasiclassical kinetic equation for the 
electron density matrix, /, which takes the form 



\ {ap(Bp • 6), 9r/p} + i[Bp • 6, /p] + eEtot • dj^^ = 



(7) 



In the above equation Ejot is the total electric field in- 
cluding both external and induced potential contribu- 
tions, and is the collision integral. [5] We allow for 
an imperfect quantum Hall effect by considering the case 
in which carriers are present in at least one of the bands 
due either to doping or to finite temperature. 

The distribution function can be decomposed into 
scalar and vector pieces, /p = rip -I- 6 • fp, and the vec- 
tor fp further separated into contributions parallel and 
perpendicular to Bp, /pbp and fp . (bp is a unit vector 
in the direction of Bp.) In this parameterization of the 



density matrix rip and /p specify valence and conduction 
band occupation numbers and fp interband coherence. 
The kinetic equation for the full density matrix can be 
separated into a set of equations for these components. 

The model's intraband response is entirely stan- 
dard [TO], except that scattering on the Fermi surface 
is influenced by the inner product of the momentum- 
dependent conduction band states. For the conduction 
band wc find that 
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dt 
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dnpiBp) 
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,<5(5p-i?pO(l + bpbpO(/^-/^,), 



(8) 



where Vp = w^p/Bp is the band velocity appropriate 
for the conduction band of Hamiltonian ([6|. It fol- 
lows that the longitudinal conductivity, is related 
to the diffusion coefficient via the usual Einstein relation 
<yxx = vpG^Dp, where vp = B-p^jl-Kv^ is the density of 
states at the Fermi level. The absence of a longitudinal 
current in equilibrium, assumed in the macroscopic the- 
ory, then follows from the cancelation between the second 
(diffusion) and third (drift) terms of the left-hand-side of 
Eq. (HI, when j<' is replaced by its equilibrium Fermi 
function value. The diffusion coefficient Dp = VpTtr/2 
with 



— 1 'f^irnp 

'''tr — 



^„2p| + A2 ■ 



(9) 



Hall response - We have seen above that even in the pres- 
ence of screening there is a residual radially symmetric 
electric potential at the surface for i — >■ oo. The purpose 
of the following calculation is to determine whether or not 
that potential can drive an azimuthal Hall current which 
contributes to the orbital magnetization. The naive guess 
that one just has to multiply the screened electric field 
with the intrinsic Hall conductivity to find the current 
fails because gradients in the density of carriers, all of 
which generally carry intrinsic magnetic moments [Sllllj. 
also yield an azimuthal current. The additional contri- 
bution can cancel the azimuthal electric field response 
either completely, as it does in the longitudinal case, or 
partially. 

Since the response we seek to evaluate includes the 
time-reveral-symmetry broken system's anomalous Hall 
effect, we should include side-jump and skew scattering 
contributions [H] to describe it fully in the presence of 
impurities. Since these are dependent on impurity scat- 
tering at the Fermi surface, they can be obtained by con- 
sidering the leading quasiclassical corrections to Eq. 
In the case of a uniform electric field, the quasiclassi- 
cal kinetic equation for conduction band electrons has 
been derived in Refs. [131 E]. This equation general- 
izes Eq. ([8| to include an anomalous distribution genera- 
tion term coming from the collision integral, and beyond- 
Born-approximation skew scattering amplitudes. Since 
we are interested here in response to a non-unform static 
electric field, we need to generalize the quasiclassical 
Boltzmann equation of Refs. [T31 Ej to the non-uniform 
case by adding a drift term, Vpdrf^, just like the one in 
Eq. Q , but with Vp now including not only the band ve- 
locity, but also anomalous and side-jump corrections. It 



is then a simple matter to see that all electric-field drive 
terms vanish in that equation in local equilibrium. There- 
fore, side-jump and skew scattering contributions need 
not be considered and the entire Hall response comes 
from the intrinsic contribution. 

The intrinsic contribution should be obtained from the 
equation for fp . Importantly, since we need not consider 
the side jump contribution, we can simply drop the con- 
tribution to the collision integral for fp coming from /p, 
since the latter gives a contribution to side-jump pro- 
cesses only. [TS] Further, for a sufficiently clean surface, 
such that BpTtr ^ 1, we can also neglect the collisional 
relaxation of fp as compared to the precession term, com- 
ing from the commutator on the left hand side of Eq. ([7| . 
The general expression for the static limit of fp is thus 
obtained simply by isolating the inter-band terms on the 
left hand side Eq. Q. We obtain 

2Bpf^ = ((Vnpap)Bp) X bp + {{e-Etotdp)ip) x bp. (10) 

The second term on the right hand side of the above 
equation leads to the standard intrinsic contribution to 
the Hall conductivity due to the interband coherence cre- 
ated by the electric field. The first term on the right 



hand side of Eq. ( 10 1 is the response to the equilib- 
rium density inhomogeneities. Its contribution to the 
current can be seen to equal the curl of the internal 
magnetic moment [5] density of quasiparticles, which is 
nonuniform in space. The right hand side of Eq. ( 10 ) 



does not necessarily vanish when equilibrium values are 
used for Up - (f^p + f^p)/2 and fp = [f^p - fSp)/2 
(/o " — np{±Bp — /i + e(j)tot))- This property contrasts 
with Eq. (|8|, in which both left and right hand sides 
vanish under the local equilibrium ansatz. 

Substituting the equilibrium values gives Vrip = 
-eEfot(dnF(Bp - pt)/dBp - dnp{~Bp - fi)/dBp)/2 and 
fp = hp{np(Bp — fi) — np{—Bp — fJ.))/2. Substituting 



these expressions in Eq. ( 10 1 and taking the local direc- 



tion of the electric field to be along the x axis we obtain: 



fp = - \eEtotBp{dp^hp)xhp ^ 



dB^ 



np{i'Bp - /i) 
(11) 

This result will recover the usual intrinsic anomalous Hall 
conductivity when the derivative acts on the ^ factor 
only. When the derivative acts on both factors we obtain 
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d ( np(yB-p — /i) 



(12) 
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FIG. 1; (Color online) The dependence of the ratio of the 
azimuthal current, je, to the radial electric field, Er, on tem- 
perature for different values of the chemical potential. 



Note that the equilibrium value of jy/E^ is not the 
Hall conductivity. The ratio instead describes equilib- 
rium currents that flow along equipotential lines of the 
screened external potential and generate a contribution 
to the orbital magnetization. 

The right-hand-side of this expression vanishes for the 
2D massive Dirac equation model for temperature T — > 0. 



In this special case Eq. ( 12 ) reduces to 



Jy_ 
Ea- 



An 



{npi-A- ^i)-nF{A- ^i)), (13) 



which vanishes for any fj. > |A| since the expression in 
brackets on the right hand side is the Fermi factor dif- 
ference between the top of the valence band and bottom 
of the conduction band. Perfect cancelation occurs be- 
tween the homogenous system anomalous Hall response 
and the current due to the curl of the internal quasiparti- 
cle magnetization density. The same cancellation occurs 
for generalized Dirac models Eq. (|6| with |p|-dependent 
velocities and constant A, as long as the p-integrals are 
convergent. This precise cancelation is however depen- 
dent on our neglect of collisional relaxation in the equa- 
tion for fp , which would lead to 1/ Artr corrections. The 
cancelation is also imperfect at finite temperature; sub- 
stantial current signal can be recovered, as illustrated 
in Fig. [T] The azimuthal current vanishes not only for 
T — > but also for T — > oo and is therefore a non- 
monotonic function of temperature. 
Discussion- When an external charge is placed near 



the surface of an ideal TI with weak time-reversal sym- 
metry breaking, it induces an azimuthal current that 
produces [S] the same magnetic field as would be produced 
by a magnetic monopole located below the TI surface. 
The axion electrodynamics model[31 13] of TI magneto- 
electric and magneto-optical properties |71 [T3J [TB] ele- 
gantly captures this intriguing property. In this paper we 
have examined how the Hall current response is altered 
by the samples imperfections which always results in a 
finite longitudinal conductivity axx = g{^^ /h). Systems 
with a finite a^x are not fully described [71 [T7] by the axion 
electrodynamics model so we develop our theory directly 
in terms of surface state electronic properties. We find 
that the apparent monopole position moves away from 
the TI surface with a velocity vm = acg. Since graphene 
based two-dimensional electron systems, which are sim- 
ilar to TI surface states, canjlj^ have g values ~ 10""^ 
or smaller when time-reversal symmetry is broken by an 
external magnetic field, there is a reasonable hope that 
it will be possible to obtain TI samples in which vm is 
small enough to enable observations in which axx plays 
no role and the axion electrodynamics model is applica- 
ble. There is a considerable recent experimental effort in 
this direction. [20] 

In the long-time limit after the external charge screen- 
ing process has been completed, we find that the az- 
imuthal current response has two contributions, one pro- 
portional to the Hall conductivity and treated previ- 
ously by Zang and Nagaosa,PT] and one proportional 
to an external potential induced change in the internal 
magnetization[TT] of the surface states. For the particu- 
lar case of a massive Dirac model the two contributions 
cancel exactly in the clean T — Q limit in the presence of 
a Fermi surface. We obtain this result using a quasiclassi- 
cal kinetic equation approach, which may not be reliable 
near band edges due to both quantum and non-linear 
screening effects, but nevertheless starkly demonstrates 
the distinction between azimuthal current and Hall con- 
ductivity responses. The special properties of the massive 
Dirac model are related to its well known [22] unusual 
orbital magnetization properties in the uniform system 
limit. In general the magnetic fiux induced by an elec- 
tron charge near a time-reversal symmetry broken TI sur- 
face is dependent on the |p|-dependence of the exchange 
potential Ap and disorder effects, and not simply on sur- 
face's Hall conductivity. 
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